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Abstract 

We reanalyze the non-linear population dynamics of a Bose-Einstein Conden- 
sate (BEC) in a double well trap considering a semiclassical approach based on a 
time dependent variational principle applied to coherent states associated to SU(2) 
group. Employing a two-mode local approximation and hard sphere type interac- 
tion, we show in the Schwinger's pseudo-spin language the occurrence of a fixed 
point bifurcation that originates a separatrbc of motion on a sphere. This sepa- 
ratrix corresponds to the borderline between two dynamical regimes of Josephson 
oscillations and mesoscopic self-trapping. We also consider the effects of interaction 
between particles in different wells, known as cross collisions. Such terms are usually 
neglected for traps sufficiently far apart, but recently it has been shown that they 
contribute to the effective tunneling constant with a factor growing linearly with 
the particle number. This effect changes considerably the effective tunneling of the 
system for sufficiently large number of trapped atoms, in perfect accord with exper- 
imental data. Finally, we identify analytically the transition parameter associated 
to the bifurcation in the generalized phase space of the model with cross-collision 
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terms, and show how the dynamical regime depends on the initial conditions of the 
system and the coUisional parameters values. 

Key words: Bose-Einstein Condensation, Two mode approximation. Non-linear 
dynamical transition. Self-trapping 
PACS: 03.75.Lm, 03.75.Kk, 03.65.Sq 



1 Introduction 

Understanding of coherent quantum tunneling of matter waves is an impor- 
tant issue in physics. In particular, the dynamics of interacting fields can 
be extremely complex, giving birth to very counter intuitive phenomena like 
coherent non-spreading wave packets and nonlinear self-trapping. These phe- 
nomena have been clearly observed in the dynamics of atomic Bose-Einstein 
condensates (BEC) in double- well and periodic potentials pl2ll3] . better than 
in superconducting Josephson junction arrays, thanks to the possibility to 
design well separated traps and due to the small dissipation in atom optical 
contexts. The wave tunneling in such system is analogous to the Josephson 
effect and results in oscillatory exchange of the condensed neutral atoms 
between the adjacent traps, as has been suggested by several authors |i] . 
The nonlinear Josephson oscillations (JO) including the many-body interac- 
tion inside each trap for the atomic BEC in a double-well trap potential are 
discussed in [GfTfS] . where the suppression of the oscillation, the so called 
self-trapping phenomenon, happens when the initial population imbalance is 
above a threshold value. Based on the modeling of the Gross-Pitaevski equa- 
tion using the mean-field factorization and two-mode approximation, Milburn 
et al [B] presented an analytical solution for the dynamical transition to the 
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mesoscopic self-trapping (MST) regime (as the number of particles exceeds a 
threshold value) in a particular initial condition where all the atoms are in one 
of the wells. More recently the possibility that the cross-collisional interaction 
become important for sufficiently large number of particles has been discussed 
in [H], hence altering the regime of transition from JO to MST. 

The purpose of this paper is first to present the complete mean field dynamics 
of the atomic BEC in two wells, complementing the work done in Milburn et al 
in] by representing this integrable dynamics in terms of appropriate coherent 
states, which leads naturally to a representation on a sphere p^llfl2] . We 
show all the relevant fixed points and its stability and discussing the physi- 
cally permitted regimes compatible with the approximations taken to obtain 
the model. This allows us to understand all the transitions between different 
dynamical regimes as a function of the collision parameters, the tunneling 
coefficient and the number of atoms in the condensate. We also find the bi- 
furcation condition for the dynamical fixed point of the mean field equations 
that determines the appearance or absence of the self-trapping regime in the 
system. Secondly, we generalize this treatment by including the collision terms 
between atoms in different wells, owing to the overlap of the tails of the con- 
densate wave functions in each trap, which is shown in [9] that, depending on 
the number of atoms, this term is non-negligible. We also show the effect of 
the cross coUisional terms on the dynamical transitions, thus determining the 
physical regime where it can be observed, a feature not realized in ref . [1^ . 

The paper is organized as follows: in Section II we quickly show the two mode 
approximation model for the Bose-Einstein condensate in a double well [6], 
but including the cross-collisional terms [9J and restate its Hamiltonian in 
terms of the Schwinger quasi-spin operators. Section III is devoted to the 
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analysis of the dynamics of a classical analogue of this model. By studying the 
semi-classical Hamiltonian obtained based on the time dependent variational 
principle (TDVP) applied to coherent states associated to SU(2) group, we 
localize all the fixed points of the model on the spherical phase space and its 
stability. By varying one control parameter of the model we show a bifurca- 
tion of one of the fixed points on a sphere, which causes the appearance of a 
separatrix of motion, dividing the phase space in two regions of different dy- 
namical regimes corresponding to the Josephson oscillations and mesoscopic 
self-trapped motions of the condensate. In section IV we show the connection 
of the classical mean field bifurcation with the quantum dynamical behavior, 
illustrated by Husimi distributions, and the spectra underlying such behavior. 
Also, we show how such a dynamical structure on the sphere is modified when 
the cross-coUisional terms are included. Finally in section V we present our 
conclusions. 

2 The two-mode local approximation (TMLA) 

Considering a trapping potential V{r) with two equivalent global minima at 
f± = ±qox, for simplicity, we assume the potential to be harmonic with fre- 
quency cu in both axis y and z and, without loss of generality, having null value 
at the minima. In this way we obtain the following single particle Hamiltonian 
with mass m: 

where b is an adjustable parameter associated to the height of the potential 

between the two wells, which we choose as 6 = ^ in such a way that the 
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harmonic approximation V±^^ = ^^y^ [{x =F go)^ + Z/^ + -2^] of the potential over 



each minimum becomes isotropic. If the distance 2go between the minima is 
large enough we can suppose that the harmonic approximation around each 
minima acts as independent trapping potential to the more populated low 
energy states of the BEC. Thus, we can suppose that the two practically 
degenerated ground states of the global potential are the symmetric and anti- 
symmetric combinations of the approximated harmonic ground states, given 



by simple Gaussians {r\u±) = 



2d2 



with width d 



as m 






X 








y 


^ \ / 

X . 

MrJ> 





Fig. 1. Left side: scheme of the double well potential with localized wave functions 
u±{r). Right side: Bloch sphere with the definitions of the angles parametrizing the 
atomic coherent states. Notice that the angle 9 is different from the one of usual 
spherical coordinates. 



Now consider a system of interacting bosons, where the Hamiltonian in a 
second quantized form [13j can be written as 

+ d^rd^r'tlj\r)ij\r)V'{f,r)ij{r)ij{r). (2) 

Then, supposing that the single level energy spectra is not considerably af- 
fected by the interaction, a TMLA is done: ^/^(r) ^ {f\u+)d^ + {f\u^)d^, where 
d± are boson annihilation operators on the states \u±). We also approximate 
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the interaction potential by a hard sphere type potential for low energy parti- 
cles : V^'(r*, f^) = Vo5{f— r') = ^^^^5(r — r'), where a is the s-wave scattering 
length [ll]. By inserting the above approximations in Eq.(|2]) and introducing 
the overlap e = = e 5^, which is small (e <^ 1), we define the tunnel- 

ing parameter Q = 2(n_|if and the collision parameters k = ^ J d^ruf{r), 
T] = ^ J d^ruf{f)u^^{f) = Ke^ and A = ^ / d^ruf{r)uj{r) = ne'l'^ for % ^ j 
= +,— ). Discarding constant terms we get the two-mode Hamiltonian 
which describes the BEC in a double well. However, to explore the natural 
group structure of the model (because of the fixed number of particles A^) it 
is more convenient to adopt the pseudo-spin operators approach introduced 
by Schwinger [15j| by defining the following operators: 

dLd_ — d\d+ (iL(i+ — d\d_ d\d^ + dLd^ 
Jx = ^ , Jy = t ^ and J^ = ,(3) 

with J = N/2. In the many-body theory Jx can be interpreted as proportional 
to the position of the condensate in the x-axis; Jy as the linear momentum of 
the condensate along the same axis and as the population imbalance be- 
tween the symmetric and antisymmetric energy eigenstates. The Hamiltonian 
(|2| in terms of these operators assumes the following form: 



H 



2A{N - 1) + ^ 



J, + 2(/t-r/)j2 + 4r7j2. (4) 



This is a Lipkin-Meshkov-Glick (LMG) type Hamiltonian, which has been 
widely discussed in the literature [T6|T7] . Here, k is the self-collision parameter 
of the condensate and, as can be seen from its definition, it is much larger 
than the so called cross collision terms rj, A (non-independent), since e <^ 1 
[6]. However, we cannot neglect these lower order parameters, as discussed in 
[0], due to its presence in the tunneling term of the Hamiltonian though the 
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effective tunneling parameter ^l' = 2[2A(N — 1) + 17/2], wliere A (altliougli 
small) is multiplied by an extra factor proportional to A^, which in a typical 
experiment can range from 10^ to 10^°. Indeed this dependence of the tunneling 
rate on the number of atoms has shown to be quite relevant for the indirect 
relative phase inference of one of the two-mode condensates through the atomic 
homodyne detection as proposed in Ref. [18]. 



3 A Classical Analogue 



In this section we obtain the semiclassical Hamiltonian of the model, based 
on a TDVP [19J via the Q-representation of the Hamiltonian ^ {9, (f)\H\9, 0) 
in the so called atomic {SU{2)) coherent states [20]. This is given by: 

where in the last equality r = e~*'^tan | with 9, (j) as defined in Fig.Q. Using 
the above expansion in the angular momentum basis, we obtain 

{9, m9, 0) ^ n{r, r*) = -J^'\^ + - v)^^^,ir* + rf 



J{2J-l)r*r 
1 + r*r r 



Then, we get the following canonical Hamilton's equation of motion for the 
generahzed variables {q,p) (with r = g+^P ^) in the phase space 

i/4J-g2_p2 



^ = O'p _ _ ^) ^^-^q'p + r^ ^'^-\ ^\ {2q^ + 2p^ - 4 J), 

-r/ ^^'^~^^ g(2g' + 2p'-4J). (7) 
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By imposing g = p = in ([T]), and transforming the canonical variables to the 
Bloch sphere angles 6*, 0, we obtain the following four sets of fixed points: 



(a) 9 = 0, (p undetermined (south pole); 



(h) 6 = 2 arctan 



R^{k -n) + 



(c) 6* = 2 arctan 



0, for g > 
7r, for g < 

7r/2, for j9 < 
-7r/2, for p > 



(d) 9 = 71, (p= < 



+ arctan ( < / ^ ^ | , for gp < 



arctan ( < / ^ '^^ ^ | , for gp > 

i2J-l „„J ^ — ^2J-1 



where we defined k = {k — 1])^^^ and n = riHrr ^"^^ ^ = ^ ^ ^ 



Now, analyzing the stability [21] of the above points and considering the pa- 
rameter regimes of the system compatible with the approximations discussed 



in Sec. II, i. e.. 



> A,r7 (A > r/) (8) 

we show that: (i) The trivial fixed point [S^a) at the south pole (9 = 0), is 
stable in all the domain of parameters allowed by condition (Isj). 



(a) Fixed points sj^b) are only defined if R'^{k — n)>^ and are stable except 



along the critical line R^{k — n) 



2 • 



(Hi) Fixed points 8j^c) exists only if R n > and are stable except if p = 
(reducing to the trivial solution). However, this is physically non-achievable 
for the above conditions (Isl). 
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(iv) All the fixed points |8](d) correspond to a single point on the Bloch sphere 
= 71 (the north pole of the sphere) and are unstable for R^{k — n) > Thus, 
when the two fixed points |8|^b) become stable and depart from the north pole, 
the fixed point becomes unstable, thus characterizing a bifurcation. 

The bifurcation condition obtained above, R^{k — n) = as expressed in 
terms of the parameters of the Hamiltonian is given by: 

(«,_3^)(Ar-i)=2A(iV-l) + ^. (9) 



4 Dynamical Transition between JO and MST regimes 



Now we discuss the quantum dynamics associated to the mean field classical 
analogue of the previous section. From the stability condition ^ obtained 
above, in the absence of cross-collision terms the bifurcation condition would 
give 

k{N-1) = ^. (10) 

This corresponds to the transition between different dynamical regimes of the 
model described in pi22] p~} In the JO regime the position of the condensate, 
given by (Jx), oscillates in time around zero, indicating no preferential tunnel- 
ing to any of the wells. In the MST regime, (Jx) oscillates around a non-zero 
value, showing that part of the condensate is mesoscopically trapped in one 
of the wells, as shown in Fig|2] 

^ Note that the transition condition is not the same as in [6] , since there it is chosen 
to fit a particular set of initial conditions, whereas here we choose to consider the 
exact point of transition of the system. 
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Fig. 2. Semiclassical time evolution of {^){9,(j)) in the MST and JO regimes for 
initial condition 9 = ^ and = 0, corresponding to all atoms initially in one of 
the wells. The parameter values are chosen N = 100, «; = ^ and the plots are in 
terms of adimensional parameter 0,t, with cross collision parameter (a) rj = k/100; 
(b) rj = k/10. 
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Fig. 3. Quantum mechanical time evolution of {6,(j)\^\9,(f)) in the MST and JO 
regimes for the same initial condition and parameters of Fig.Q as a function of fit. 
(a) T] = k/IOO; (b) rj = k/10. 

In Fig.([2| (a) and (b) we show the semiclassical time evolution of {^){t) for 
the same parameter values corresponding to the quantum mechanical results 
in Fig.(|3])(a) and (b) respectively. It is clear that for r] = k/100 the initial 
state is related to a MST orbit whereas for = k/10 it is related to a JO 
orbit. The semiclassical evolution is very similar to the quantum one, except 
for a modulation in the quantum case which produces a sequence of collapses 
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and revivals (for r] = k/10, see the inset of Fig.(|3])(b)). 



We can better visualize the orbits of the system in the unit sphere whose 
surface is parametrized by the polar angles 9 and 0, the so called Bloch sphere. 
Considering > and k > 0, we have a fixed point at 6 = tt (north pole) 
which is stable whenever {K—3ri){N—l) < 2A(A^— 1) + ^. When this inequality 



is reversed, the north pole becomes unstable and two new fixed points appear 




Fig. 4. Trajectories on the Bloch sphere for several values of self-collision, and 
without cross collision term (r/ = 0). Fixing N = 100 and 0, = 1, the self-collision 
parameters of each plot is given by: (a) k = 0, (b) k = (c) k = ^ and (d) 



Notice that in the case rj = 0, for increasing values of k, before the transition 
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point Kc = 2{N-i) there are solely orbits associated to the JO of the condensate. 
When the self-collision reaches the critical value, the bifurcation happens and 
orbits of MST orbits appears, occupying an increasing area in the phase space 
of the system as k, increases. As already noticed in previous publications for 
particular initial conditions [6] , the increase of collision rates between bosons in 
the same site are responsible for the phenomenon of suppression of tunneling. 

The Q-representation of a Hamiltonian gives the exact classical limit of a sys- 
tem in the N —>■ oo limit [23], thus the bifurcation condition on the parameters 
{k — 37]) (N — 1) = 2A(A^ — 1) + ^ represents in this limit the exact transition 
condition between the two dynamical regimes of the model. 



Now, we consider the effect of the cross coUisional terms on the population 
dynamics, and we find a non-negligible change in the effective tunneling even 
for a number of particles not so large (A^ = 100) as shown in Fig.(|5]). 




Fig. 5. Trajectories on the Bloch sphere, considering cross collision terms with in- 
creasing values of rj and fixed parameter values N = 100, 17 = 1 and k = j^. (a) 

In opposition to the self-collision term, the increase of rj suppresses the pres- 
ence of the MST orbits from the phase space. Even for values of i] much smaller 
than K we observe a significant change in the dynamical behavior of the sys- 
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tern; namely, for solely 100 particles it is already noticeable the modification 
in the effective tunneling of the system. This fact goes in the same direction of 
the experimentally observed situatiorp] p] . Notice that in this work we always 
keep the cross-collision rates considerably lower than the self-collision rates, 
respecting the hypothesis used to obtain the model, and we select only the 
physically allowed bifurcations from those listed in [H] obtained from Eq.Q. 

In order to show the consistency of the quantum dynamics with both dynami- 
cal regimes, we show in Figs. ^ and ([T]) a sequence of plots of Husimi functions 
on the sphere for time evolved states with the same initial coherent state, but 
for parameter values corresponding to MST and JO regimes, respectively. 

Although the initial state for Figs ^ and ([T]) is the same coherent state on 
the sphere centered at the position where the x-axis crosses the unit sphere 
(represented in Fig. ^ for ^It = 0) the time evolved states are completely 
different due to their respective dynamical regimes. Fig. ([6]) corresponds to 
the MST regime as shown in Fig.(|3^) with the classical trajectories in phase 
space corresponding to the one shown in Fig([5^). Thus, the time evolution 
shows first a spreading of the Q-function along the classical trajectory {Qt = 
5, 10) in the trapped region of the phase space, and then shows effects of self- 
interferences {Qt = 20, 30) forming mesoscopic superpositions with three and 
two well defined peaks, and finally re-constructs the one-peak quasi-coherent 
state at Qt = 63. Notice that Qt = 30 corresponds to the time at which a small 
recovering of oscillations (revival) happens in the middle of the collapse and 
= 63 to the revival time in Fig.(|3^). This behavior is a clear demonstration 

^ This does not exclude other possible explanations for these effects such as the 
one presented in [25] . 
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Fig. 6. Plot of a time evolved Husimi functions for MST initial condition with 
parameter values N = 100, rj = k/100 and k = 1/N for nt = 0, 5, 10, 20, 30, 63. 

of the tj^ical collapse and revival of the condensate state phase. 

In Fig.Q on the other hand, it is shown the evolution corresponding to the 
JO regime as shown in Fig.([3|D) and the classical trajectories in phase space 
corresponding to the one shown in Fig([5]3). The initial distribution at f2t = is 
the same as in fig. Q. One can see the spreading (collapse) of the distribution 
surrounding all the sphere characterizing the JO {Qt = 5, 10), and then the 
self-interferences also happening but in a less coherent way {Qt = 30,60). 
Finally at Qt = 250 the best re coherence of this case is shown, but clearly 
in a much less coherent reconstruction, as expected due to the much longer 
classical trajectory to be traversed, and this is consistent with the less perfect 
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Fig. 7. Plot of a time evolved Husimi functions for JO initial condition with param- 
eter values N = 100, t] = k/10 and k = 1 /N for nt = 5, 10, 30, 60, 250. 
revival in the inset of Fig. (pi; 



Another signal of the dynamical phase transitions (JO — > MST) can be seen in 
the energy spectra of the Hamiltonian (|4]). In the spectra shown in Fig.([8])(a) 
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Fig. 8. Energy spectra En as a function of the label n for fixed parameter values 
N = 100 and $7 = 1. (a) without cross-collision terms (77 = 0), and several values of 
self-collision parameter k. (b) with cross- collision terms, fixing 1^ = jj and several 
values of rj. 



we observe that for self-collision parameter above the critical value there is 
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the appearance of an inflection point ([23]) and a doublet structure. The term 
responsible for the doublet structure is the term in the Hamiltonian Q (as 
shown in the inset) representing the self-trapped states when it is dominant 
(namely, eigenstates of well deflned difference of the number of particles be- 
tween the two wells). In Fig.(j8])(b) we show the spectra with flxed = jq and 
several values of cross collision parameter rj. As we increase rj (but keeping k 
significantly larger than r^), the inflection point and the doublet structure are 
clearly suppressed. 

5 Conclusion 

To conclude, we confirm first that the effect of cross-collision in the two mode 
approximation model for the dynamics of the Bose-Einstein condensate in a 
double well trapping potential should not be neglected in all circumstances 
|9]. This is so because the intensity of this effect on the effective tunneling 
of the system increases linearly with the number of trapped particles, which 
in a typical experiment can easily reach a mesoscopic amount of atoms. Such 
effects has been shown in the simulations presented in Figures 
these are qualitatively in accordance with existing experimental results [2]. 

We also showed that the transition of the dynamical regime of the system is 
associated with a bifurcation in the generalized phase space in a semiclassical 
approach of the problem, in such a way that the dynamics of the condensate 
- Josephson oscillation or macroscopic self-trapping - is extremely dependent 
on the initial conditions chosen for the populations of the wells and its phase 
relations. Our semiclassical approach using PVDT and SU(2) coherent state 



4||7p) and 
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guarantees both the correct phase space topology and the particle number 
conservation. Also, taking advantage of the SU(2) group structure of the model 
leads naturally to a phase space motion on a sphere, facilitating the analy- 
sis and visualization of the dynamical structure associated to the nonlinear 
quantum dynamics of the double well condensate, namely JO or MST. 

Finally, we show a situation where the cross-collision terms can even suppress 
the bifurcation and thus the appearance of the self-trapping regime in the 
Bose-Einstein condensate in a double well. 
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